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Abstract 

A new method is introduced for probing Topological Superconductors. 

The integration of the superconding fermions generates a topological Chern — Simons sound 
action . 

Dislocations induce Majorana zero modes inside the sample, resulting in a new Hamiltonian 
which couple the Majorana modes, the electron field and the non-Abelian strain field sound. This 
Hamiltonian is used to compute the anomalous sound absorption. The Topological superconductor 
absorbs sound below the superconducting gap due to the transition between the quasi particles 
and the Majorana fermions. 

The sound waves offers a new tool for detecting Majorana fermions. 
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I-Introduction 


Topological materials have been discovered recently BE]. A new class of insulators 
coined Topological Insulators are characterized by the second Chern number [TB] which 
can be measured using the Faraday and Kerr rotation [2j. Similarly to Topological Insu¬ 
lators, Topological Superconductors have been found and identihed by the presence of the 
Majorana zero modes [IHZl E]. An indirect observation for the Majorana fermions was 
obtained from the differential conductance measurement [8]. Since the electrical current is 
not conserved in a superconductor, the electromagnetic field can not be used to identify the 
Topological Superconductors. The response of a Topological Superconductor is character¬ 
ized by a topological Chern — Simons sound action as we have for the Yang-Mills action 
[3]. The coeficient of this action is the central charge c which counts the number of the Ma¬ 
jorana edge which are the signature for the Topological Superconductor. The topological 
action contains high derivatives and it is unlikely to be obseved in the laboratory. For this 
reason we need to hnd an alternative method. 

Sound waves have played an important role in the field of Superconductivity. The first 
experimental measurement of the superconducting gap was by ultrasound attenuation eg. 
A typical ultrasonic attenuation wave has a frequency of 10 — lOOMifZ which is too small 
to break a pair and excite a quasi-particle, since gaps are in the range of milli-electron 
volts. This suggests that in order to enhance the sound absorption one needs to increase 
the density of the normal matter in the superconductor. This can be done with the help 
of vortices or dislocations. 

In a superfluid it has been shown that in the mixture ^He-^He, the density of the normal 
matter pn increases with ^He [niEo]. The technique to measure the normal density pn is 
based on the torsional oscillator [18], the period of the oscillation is related to the moment 
of inertia of the normal density. Similar measurements have been done by HU to investigate 
^H,-B . 

Here we would like to propose a new method based on the coupling of sound waves 
to Majorana fermions and electrons. We will present a theory which demonstrates that 
the presence of the induced Majorana modes affects the response function . We compute 
the sound polarization and obtain the normalized sound wave equation which is used to 
evaluate the shear and longitudinal stress. The stress serves as a probe for the Majorana 
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fermions. The stress field generates dislocations and strain fields. The effect of the stress 
held on the superconductor can be measured using impedance techniques, by measuring the 
change of the quality factor Q and resonance frequency of an ac-cut quartz transducer which 
oscillates in a shear or longitudinal mode 

Topological Superconductors are characterized by a winding number and the central 
charge c in the the region where the chemical potential is positive. At the interface between 
the region of positive chemical potential and a negative one, (where the superconductor is 
not a topological superconductor) zero modes (half vortices) appear [T2HI1]. Half vortices 
are identihed with the Majorana zero modes. The formation of the Majorana fermion is due 
to the vanishing of the Topological Superconducting density, which on a two dimensional 
surface is identihed with the vortices. Solids are characterized by dislocations which are 
either induced by large external deformations or are a result of the crystal growing process. 
On a surface, the presence of the dislocations are similar to two dimensional vortices. The 
Majorana zero modes are bound to the vortices on the surface of a Topological Supercon¬ 
ductor. This suggests that dislocations control the density of the Majorana fermions. By 
increasing the density of the vortices in a Superconductor, one can enhance the acoustic 
absorption. We hnd an anomalous absorption for frequencies A — ea < oo < A + ea, A 
is the superconducting gap and €„ represents the overlapping energies for two Majorana 
fermions. The anomalous absorption occurs in the forbidden superconducting frequency 
region u < 2A. The anomalous absorption represents the hnger print of the Majorana 
fermions in a Topological Superconductor. 

This results are obtained with help of a new Hamiltonian. Using a space covariant 
transformation induced by the sound waves we derive the coupling between the non-Abelian 
strain held, the Majorana fermion and the electron held. The sound waves ’’deforms” the 
Topological Superconductor. The deformation is obtained with the help of the coordinates 
transformation method [2IH23]- 

The plan of this paper is as follows . In Sec. II we introduce the Topological Super¬ 
conductor Hamiltonian. Using the p — wave model we construct the deformation caused 
by elasticity, following the methods given in the literature [I51I2IH23]. We show that the 
integration of the fermions generate the topological Chern — Simons term with the central 
charge which counts the Majorana modes on the boundary of the sample. In the second 
step we consider dislocations which induce Majorana zero modes inside the sample. For 
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this situation we obtain a new Hamiltonian and show that the Majorana modes couple to 
the electron held and sound. The details are given in Appendix -A and Appendix -B. 
Sec. Ill is devoted to the computation of the sound absorption, the details are given in 
Appendix -C. We compute the shear and longitudinal stress and identify the anomalous 
sound absorption which conhrms the presence of the Majorana modes. In section IV we 
present our conclusions. 

II- The Topological Superconductor in the presence of an elastic strain field 


In order to demonstrate the anomalous absorption we consider a p — wave superconductor 
HE]. The p — wave superconductor is engineered using materials with strong spin orbit 
interaction and magnetic Zeeman heldsi which are in proximity to an s-wave superconduc¬ 
tors [I21E3] . In order to study the stress response and absorption we need to know the 
coupling between the sound waves and the superconductor. The form of the electron phonon 
Hamiltonian is problematic for spinors. The symmetry of the coupling is obtained from the 
coordinate transformation method [T3] . When a sound waves excites the crystal, the crystal 
coordinates are modihed. The external strain held modihes the position of the crystal atom. 
As a result the p — wave is replaced by a deformed superconductor Hamiltonian: 

H = _|_ ^jj{cr-p-wave) _|_ jj(cr) 


fjip-wave) 

is the model for the p — wave superconductor in the absence of the sound waves. 


^{p wave) ^ ^ j d‘^rC\f,t) -iT'^'^{di+id2) 


( 2 ) 


A(r, f) the pairing order held, prir) is the space dependent chemical potential and r^,r^ 
are the Pauli matrices in the particle-hole space. We assume that in one region Pf{I) > 0 and 
in the complimentary region < 0. For Pf{I^ > 0 the superconductor is topological 

and is characterized by the topological invariant with the Chern number QgZ [T2], in 
the region Pf{I) < 0 the superconductor is non topological. At the interface Pf{I) = 0 
(between the two regions) the spectrum will contain bound states, Majorana zero modes. 
The change of sign of the chemical potential in space gives rise to the Majorana fermions. 
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In the literature a few proposal exist, for the Majorana fermions, mostly at the interface 
with a Ferromagnetic region [T2] and vortices HU- Recently new proposals based on the 
interplay of magnetism and superconductivity P dU] have been introduced. The change 
of sign of the chemical potential in a p — wave superconductor can be achieved by the 
presence of vortices or dislocations (introduced by the growing process of the crystal) or 
large mechanical stresses. The presence of the dislocations localized at position Rb gives rise 
to discrete points where the chemical potential ^ Rb) vanishes (see Appendix-A). 

The crystal Hamiltonian controls the sound propagation. The crystal Hamiltonian 
is given by: 


^(cr) _ / ^2^ 


-^2 

'— + - ( diu^ + diu^ 
2p 2 


2 A 


i\2 


-,hj = 1,2 


(3) 


The solid Hamiltonian contains the kinetic and potential energy. The kinetic energy 
-*2 

^ is determined by the canonical momentum n and mass density p. The potential energy 
is characterized by the Lame elastic constant A and shear modulus p m 
fextX^R) = is the force of an A — C quartz transducer which pumps energy into 

the solid. The external source Hamiltonian jg given by: 

= j cfrf^^t.irR) ■ u{r,t) (4) 


The sound propagation modihes the p — wave .In the presence of the sound waves, the 


coordinates are modihed from f=x= to x + u{x,t) = X = 


X^=\X' 


a=2 


where 


u{x, t) = t),u^‘^\x, f)]. The sound modihes the term r^'da (in the Hamiltonian Eq.(2)) 

to Hg, is computed with the help of the transformation X{x) induced 

by the sound. As a result a non-Abelian strain held A{x,t) is generated. The non-Abelian 
held A{x,t) is written in terms of the elastic strain held.The non-Abelian form is due to 
the particle-hole operators r^,r^ represents the ehect of the sound waves which is written 
in terms of the elastic held u{x,t). The full derivation of }jA(^f°rmed-p-wave) Hamiltonian is 
given in Apendix -B. The integration of the fermions with Hamiltonian 
and pf > 0 generate a topological sound action. 




—sound 


967r 


dt / d^r 


- dkUlj) + -e'"’’'’‘'uJi^aUJi^bUJi^c 


(5) 


c counts the number of the Majorana modes on the edge of the sample for pp > 0. Due to 
the high order derivatives is unlikely to be observed in the laboratory. 
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We will consider a chemical potential which changes sign inside the sample due to vortices 
or dislocations. Due to vortices or dislocations we will have additional Majorana zero modes. 
We will focus on the hrst non-trivial = J^{deformed-p-wave) _ j^{p-wave) 

which contributes to the absorption. Such a term is linear in the Majorana and in the 
electron held. 


^{deformed—p—wave) jj{p—wave) ^ ^^{cr—p—wave) / ^ 2 ^ 




A{f,t) = Ibo{f,t) + ribi{fA) + r 2 b 2 {r,t), A = |A|e*“, bo{f,t) = dtu{fA) 

bi{f,t) = '^|A|( cos(a)92'u(r, t) — sin(a)9i-u(r, t) j b 2 {f,t) = '^|A|( cos(Q;)9i'u(r, t) + sin(a)92w(D t) 


( 6 ) 


The Hamiltonian is diagonalized with the help of the Bogoliubov - deGennes 

transformation. The Hamiltonian has in addition to the positive eigenvalues also Majo- 

r ^ ^ 1 ^ 

rana zero modes mil!- The spinor for non-zero energies are given by \U{k),V{k)\ and 


the zero modes spinor are given by 


The Bogoliubov — deGennes 


Hamiltonian contains pairs of momentum [fc, —k\ and therefore the momentum integration 
is restricted to half of the Brillouin zone. To cover the entire Brillouin Zone we will replace 


U{k),V{k) by a four component spinor <h(fc) = U{k),V{k),U{—k),V{—k) 


1 T 


, similarly 


for the Majorana spinor we will replace t/fe(r, 0), 0) with a four component spinor 

m{f, 0) = k(r, 0), 0), U,{f, 0), U^(f, 4>) 


C(r) = ^ 
C^(r) = 


d‘^k 

■ d'^k 

(27r)2 


2n 


Akf 


—ikr 


g{k)^^k) + r]{-k)T^*{k)\ + 0) 

6=1 

2n 


b=l 


(7) 


C(r) and Ct(f) are the quasi particles operator for the superconductor. The matrix T 
ensures the (pseudo) reality conditions [TT] . 

rCt(r) = C(r), T^l{k) = ^_E{k) (8) 

The operators C(r) , C^(r) also contain the Majorana zero modes =176^4(^0) , 
(g 0)]^ where jb are the Majorana operators. Using the mode expansion given 
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in Eq.(7) we will replace the Hamiltonian ^ by: 


^jj{cr—p—wave) _ 

(Pk 


2n 


dr 


—ikr 


(27r) 




6=1 


+ h.c. 

(9) 


III- Computation of the anomalous sound absorption 


In order to perform the computation in Eq. (9) we need to know the spatial wave function 
for the Majorana fermions. Due to their localization in space and low density we can assume 
a model of randomly distributed Majorana fermions. We will consider a situation where the 


correlation within the nearest neighbor pair 


76=2a-l, lh'=2a 


is signihcant and negligible 


otherwise. This allows to introduce the fermion operators CliCa ; 


72a-1 — 




Cl + Ca 


72a = 




Cl-Ca 


, a = 1 . 2 ..n 


( 10 ) 


The overlap between the pair 


76=2a-l, lh'=2a 


is given by the overlapping energy e^. The 


Majorana Hamiltonian is given by = ^“~”eaClCa- The ground state with the 

Majorana fermions is given by |G, 6 ^= 1 ...ea=n) • This ground state |G, 6 ^= 1 ...ea=n) annihilates 
the operators rj{k) and Ca, rj{k)\G,ea=i...ea=n) = 0, Ca\G,ea=i...ea=n) = 0. We compute 
the expectation value with respect the ground state \G,ea=i...ea=n) and hnd the sound 

polarization 11 [^,a;]. 




^ /To, dt5H(‘=-'>[diuHr,t),dtui{r,t)] 


( 11 ) 


As a result the sound wave Hamiltonian is replaced by: H = i/hi") _|_ //(ea;t.) _|_ §fficr) 

is given by Eq.(ll). (It gives the sound absorption induced by the transition between 
the superconductor quasi-particles with energy E = hv\J{e{k) — fipY + |ApA ;2 the 
Majorana fermions Cq.) We evaluate Eq.(ll) to order ^ and obtain the polarization diagram 

n[g,ca]. 

/ 72 7 r ^ 

j ■^U[q,u]u{q,u)-u{-q,-u), H [g, w] = ReH [g, ca] + ImH [g, ca] (12) 

The polarization n[^, ca] is computed in Appendix-C. Using typical values for the sound 
waves frequency ca = lO^R;^ and Fermi momentum kp = 10®m“^ we hnd that the value of 


7 



the self energy 11 [^, u] is in the range 100 therefore the effect on the phonon frequency 

is negligible, and we will consider only the imaginary part of the polarization 11 [^, 0 ;] which 
determines the absorption. 

Il[qM = ^ [1 - 2|A|M cos[a - m] - 4|AP(^)"1 ■ E 

Stt L |a;| |a;| J “ 

r dE - Ef 

^Je^ - A^k^p {^Je^ - - Ef + 

0M ((0[^]0[eg - 0[-.E]0[-eg)5[a; - {E + ej] + (0[-E]0[eJ + 0[E]0[-eg)5[a; - {E - ej]) + 
0[-a;] ((0[^]0[q - 0[-.E]0[-e,])<5[a; + {E + e,)] + (0[E]0[-e,] + Q[E]Q[-iaMu + {E - e,)]) 

(13) 

In Eq.(12) we used the dehnitions e{k) = i{k) = ^,/lp = E = ^ = 

^{i{k)-f^p)2+\A\^k^. 

The hrst line of Eq.(13) is determined by the structure of the sound held A{r,t) given 
in Eq.(3). X)a=i represents the sum over the uncorrelated Majorana fermions Ca-The 
second line of Eq.(13) contains the density of states ^^d the coherence factor 

y i?2 —A2fc|, 

( introduced by the spinors ^{k), Wa{r,(l>)) 

The third line with 0[a;] {u > 0) describes the sound absorption. 

The first term in the the third line describes the creation of a Majorana hole with the energy 
Cq (destruction of Majorana particle with the energy —e^) and the creation of a electron with 
energy E. As a result the absorption of sound will occur for frequencies a; > E + with the 
threshold absorption u > \A\kF + ia- In the absence of the Majorana term the absorption 
will occur u > 2^,resulting in a forbidden frequency region u > 2|A|A;i7’. 

The second term in the third line describes a situation where the Majorana states are 
occupied at the energy e^, and the sound absorption will excite the electron to the state 
energy E. As a result the sound will be absorbed for frequencies cu > E — ia ,resulting in 
the threshold absorption u > \A\kF — ia- This absorption will be controlled by the Fermi 
Dirac occupation function which at hnite temperatures T replaces 0[—Ca] with nplea] = 

—. The absorption region \A\kF — ia < OJ < \A\kF + ia will be called the anomalous 
absorption. Due to the fact that |A|A;ir >> ia the presence of the Fermi Dirac function 
ni?[ea] in the absorption formula gives the possibility to determine the energy distributions 
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e[E-AkF_ ■ 



of the Majorana fermions. 

The fourth line with 0[— ta] (u < 0) describes the sound emission. 

We will use the polarization operator n[g, ta] obtained in Eq.(13) to compute the sound 
waves propagation in the crystal. We note that due to the random distribution of the Majo¬ 
rana fermions the sound polarization is isotropic and will affect the sound waves equations in 
a symmetric way (see Eq.(14)). The crystal Hamiltonian is normalized by the sound polar¬ 
ization n[g,a;]-u(g,cu) ■ u{—q, —u). The sound wave equation for the external force 
is given by: 


(—+ l^q^ + n[g', -|- (A -|- /r -|- 


U[q,u_ 


)qiqj u^^\q,uj) 


/W(g,cu); i,j = l,2 (14) 


We consider a situation where the force in the x direction is zero f^^\q,u) = 0, and the 
force in the y direction is f^‘^\q,u) ^ 0. We compute the shear stress the 

longitudinal stress 0 - 2,2 [ 9 , 1 ^], the strain tensor 61,2 [ 5 , 0 ;], 62,2 the impedance Z[g, co] 
in terms of the shear modulus y and Lame elastic constant A .Using the solutions obtained 
from Eq.(14) we End: 


c^i,2[g,w] 
0-2,2 [g,c<;] 

Zi,2[g,w] 

Z2,2[g,t^] 


ei,2[g,a;]2/i, ei,2[g,a;] = - iq2U^^\q,uj) + iqiu^'^\q,uj) 


e2,2[g,t^](3A-h 2/i), e2,2[g,t^]= iq2U^^\q,u 

c^i,2[g,w] _ _ r- 1 _ 

n( 2 )\^ -1 = 1^1,2 9)1^ + 2 g, CO — 

^2,2[g,co] _ _ _ 

= R2,2[g,coJ *X2,2[g,coJ = 






CO 


/ig^ 4- ilmn[g, co] 
f(3A -|- 2/i)g2 
/xg^ 4- ilmn[g, co 


(15) 


We plot the real part of the impedance Ri, 2 [^, co] for a fixed Majorana energy, ignoring the 
energy dispersion of the Majorana fermions. Due to the isotropic form of the polarization 
the longitudinal impedance R 2 , 2 [g, co] has the same frequency dependence as Ri, 2 [g, co] (see 
Eq.(14)). We will use in the plot a narrow energies distribution ^ ^ 0.05. 

The thin line gives the anomalous absorption for the Maiorana 

energy ^ 0.05 and the superconducting gap ~ 0.25. up = ^ is the corresponding 

frequency for the Fermi energy. 

The thick line shows the absorption in the absence of the Majorana fermions for frequencies 
^ (The form of the gap A|fc| is a result of the Hamiltonian given in Eq.(2).) Since 
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(Rq^ 12 ji j ,e/w^ ] 



LOjCDp 


FIG. 1: The real part of the impedance is R[q, w] . The thin line gives the anomalous absorption 

for the Majorana energy ^ ~ 0.05 and the superconducting gap 
PS 0.25. The anomalous absorption is plotted for low temperatures such that n^lea] ~ 1- The 
thick line shows the sound absorption in the absence of the Majorana, the absorption starts for 
frequencies — > =0.5 

disorder is present, it raises the question how to differentiate the absorption between the 
Majorana fermions and charged impurities. Charged impurities will give rise to absorption 
for frequencies u > E — iimp. and temperatures KpT > 2|A|/ci?, contrary to the absorption 
induced by the Majorana fermions at T —)■ 0. 

IV-Conclusions 

To conclude a new proposal for probing superconductors with sound waves has been intro¬ 
duced. An anomalous absorption is found in the frequency region \ A\kF—ea < oj < \A\kp+ea 
and the energy distributions of the Majorana fermions can be determined. This gives direct 
evidence for the presence of the Majorana fermions. The Majorana fermions are the huger 
print of the Topological Superconductors. Therefore by measuring the anomalous absorp¬ 
tion we have identihed a method for detecting Majorana fermions and probing the response 
of the Topological Superconductors. 
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Appendix -A 


Due to the large deformation introduced by the dislocation localized at Rb, the crystal 
coordinates are modihed from f = [x,y] to R{r) = [A(r'),D(r)].This allows to introduce 
the static strain held e“ = diX^, a = 1,2 and i = 1,2. Due to the dislocation the the 
area dXdY is modihed to J[Rb\dxdy where the J[Rb] is the Jacobian of the transformation, 
As a result the chemical potential yp becomes fip —t ypif) = ypJlR}]- For an edge 

(y) 

b _ 

(x) • 


dislocation the coordinate transformation is given hy \ X = x and Y = y — ^^tan 


5 ( 2 ) 

is the Burgers vector in the y direction. This transformation gives for the Jacobian 


J[Rb] = (elel - elej ) = 1 - 

yp{r Ki Rb) vanishes. 


x—R 


(a:) 


2 ^ 


As a result the chemical potential 


Appendix- B 


The sound waves held u{x, t) change the coordinates from r = x= [x,y] to x + u{x, t) = X = 
where u{x, t) = t),u^‘^\x, t)]. The strain held are dehned by El = 

e“ = are related. We have: ELi Ea=i ,2 = gij and Ei^b = QigEl For 

the time component we have El = e\5i^f 

The derivatives transform like vectors, da = daU^{x, t)di + daU^{x, t)d 2 - For the sound waves 
u{r, t) we have: e} = 5i^i — div}{x, t), ef = 62,1 — diV?{x, f), e“ = —dtu°‘{x, t) for a = 1, 2. The 
integration area element d‘^x in Eq.(2) is multiplied by the Jacobian J = 


1^2 


ete 




plp2 


Det[e“]. The deformed p — wave Hamiltonian p wave) form: 


^{deformed—p—wave) _ _ 

“ 2 


<fxDetle‘]CHf,t) [t=>( A(^£);bjv.V,) - ixf{x) ] - A(x,t) 

a=l 


+A* (x, t) E+ iY {E\Xi - iE\Xi) C(x, t) 


Vj is the covariant derivative given in terms of the spin connection: 

r^], Vj = r^]. The spin connection has been derived in [22] in terms of El 

and e“. The spin connection is determined from the zero torsion condition, Vje“ — = 0. 

We have. 


wf = 




Ei{d.E^,b 


djEi^b) 


2^ 


b,c,d 


{EiE)dtE,Pe1 


(17) 


Where Eij, ^ gijEl. 
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Once the spin connection is know we can perform the path integral over the Dirac’ 
fermion. As for the Yang-Mills theory the fermion integration in 2 -|-1 dimensions generate a 
non-Abelian Chern-Simons term. We perform the path integral integration tor the fermion 
held C\r, t) and 0(r, t) and obtain the effective sound action [ 23 ]. The topological 

term is given by, 




—sound 


967r 


dt / d^r 




dkCUj') Y e 


(18) 


where Ui^c = c counts the number of the edge modes. For the Topological Super¬ 

conductor we have pi? > 0 and c is non zero. For this case Majorana modes are on the 
edge of the sample. As a result the effective sound action allows to identify the 

Topological Superconductor. Due too the high order derivatives it is difficult to observe such 
a term in the laboratory. For this reason we will look for an alternative way to identify the 
Topological Superconductor. This result is similar to the gravitational Chern-Simons term 
[3]. 

In the presence of dislocations we expect to have additional Majorana zero modes. There¬ 
fore the sound waves can couple to the Majorana fermions and electron field. In order to 
study this effect we will consider small deformations |■u(r, f)| << a (a is the lattice constant). 
We will keep only first order terms in the coordinate transformation and neglect the effect 
on the metric tensor and spin connection 


j^{deformed-p-wave) ^ ^ A 

^2/ 


+A*{r,t)(P +iT^){E\d, 


C\r) {Eldi + iEldi) 

iE\di) C{r,t) 

(19) 


We dehne ^ wave)._ 


^^{cr—p—wave) _ ^{deformed—p—wave) _ ^{p—wave) 


The Hamiltonian ^ is given in Eq.(3). 


( 20 ) 


Appendix-C 

We use Wick’s theorem for the Green’s functions, G{k, t) = 
—z(G, 60 = 1 ...ea=„|T(? 7 (fc, r)? 7 l(/c, 0))|G, 6 ^= 1 ...e(j=„) for the Superconductor and 
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ga{r) = -i{G,ea=i...ea=n\T{Ca{t)C{0))\G,ea=i...ea=n) for he Majorana fermions. The 
spinors Wa{f, 0) are used to compute the ’’coherent” factors of the self energy n[g, ca]: 




Stt 


1^ „„„r„ A\ \ \ 2 /\^ \2 


n[g, cn] = 1 - 2 |A|f^ cos[a - m] - 4|A|^(fA) 


ca 




dE 


- A^kl 


E^ - A^kl - Ef 


n « 

- / dr ({G,ea=l.:ea=n\T(T]^{E,T)T]{E,0)^\G,ea=l:.ea=n) 

'E^ - A^kl-Ey + A‘^k^^ a=iJ 

+ {G, ea=l---ea=n\T (^ri(E,T)riyE, 0)^ |G, ea=l---ea=n)^ (^{G, ea=l---ea=n\T ^Cl('r)Ca(0)^ |G, ea=l---ea=n) 

+ {G, ea=l-..ea=n\T ^Ca('r)Cl(0) j ea=i...ea=n) j 

( 21 ) 


T = ti — t2, T stands for the time order product. In Eq.(17) we have used the notation 
E = f-; ia = y"] tan[/3{q)] = We integrate with respect t = ti — t 2 and obtain: 


n[ 9 ,c] = ^ [l - 2 |A|M eos[a - / 3 ®] - 4 |AqM )2 


Stt 


dE 


E^ - A'^kl 


U E^-A^kl-Ef 


n 

n- 


0 [^] 0 [e. 


+ 


0 [^] 0 [ea 


(y^E2 _ a2A;2, _ E )2 + A2A;2 a=i + (E + Q - ix u - {E + Q + ix 

0 [-^] 0 [e« 


^ e[-E]Q[-ia] ^ Q[-E]Q[-ia] ^ e[-E]Q[ea] ^ 

u — {E + ia) — ix u + {E + ia) + ix u + {E — e^) — ix u — {E — ia) + ix 

, a; —)■ 0 


0[E]0[-eJ ^ 0[E]0[-e, 


u + {E - ea) - ix uj - {E - ea) + ix 


( 22 ) 


At hnite temperatures the steps functions are replaced by the Fermi Dirac ocupation 
functions, Q[-E] = upiE] = 0[-ea] = nir[ej = ^^i^,Q[E] = 1 - 0[-F], 

0[ea] = 1 - 0[-ea]. 
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